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> 

O : We analyze PLANET collaboration data for MACHO 97-BLG-41, the 

Q . only microlensing event observed to date in which the source transits two 

. disjoint caustics. The PLANET data, consisting of 46 V-band and 325 I-band 

Q\ . observations from five southern observatories, span a period from the initial 

alert until the end of the event. Our data are incompatible with a static binary 
lens, but are well fit by a rotating binary lens of mass ratio q = 0.34 and 
angular separation d ~ 0.5 (in units of the Einstein ring radius) in which the 



c/3 . binary separation changes in size by Sd = —0.070 ± 0.009 and in orientation by 

66 = 5°.61 ± 0°.36 during the 35.17 days between the separate caustic transits. 



We use this measurement combined with other observational constraints to 
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derive the first kinematic estimate of tlie mass, distance, and period of a 
binary microlens. Tlie relative probability distributions for these parameters 
peak at a total lens mass M ~ 0.3 M© (M-dwarf binary system), lens distance 
-Dl ~ 5.5 kpc, and binary period P ~ 1.5 yr. The robustness of our model 
is demonstrated by its striking agreement with MACHO/GMAN data that 
cover several sharp features in the light curve not probed by the PLANET 
observations, and which did not enter our modeling procedure in any way. 
Available data sets thus indicate that the light curve of MACHO 97-BLG-41 can 
be modeled as a source crossing two caustics of a physically-realistic rotating 
binary so that, contrary to a recent suggestion, the additional effects of a 
postulated planetary companion to the binary lens are not required. 



Subject headings: binaries: general — gravitational lensing — planetary systems 



"Based on observations at: Canopus Observatory, Tasmania, Australia; Perth Observatory, 
Bickley, Australia; the South African Astronomical Observatory, Sutherland, South 
Africa; the European Southern Observatory, La Silla, Chile; and the Cerro Tololo 
Inter-American Observatory, Cerro Tololo, Chile. 



-3- 



1. Introduction 



One of the first microlensing events discovered was shown to be due to a binary lens 
( |Udalski et al. 1994D , not surprisingly since most of the hundreds of events discovered 



toward the Galactic bulge are believed to be due to normal stellar lenses in the disk and 
bulge of the Milky Way (cf., |Kiraga fc Paczyhski 199^ ; |Zhao, Spergel fc Rich 1995| ; |Alcock| 
et al. 1997|) and a sizable fraction of stars reside in binary systems. However, not all binary 



lenses can be recognized readily as such. For typical photometric sampling and precision, 
most light curves generated by double lenses will be indistinguishable from single-lens 
events if one lens component is sufficiently less massive than the other or if the components 
are separated by angular distances sufficiently smaller or larger than the Einstein ring 
radius. Nevertheless, several binary microlens events are recognized every year (cf. [Alcock 



et al. 19^ ), generally through sharp peaks in their light curves that betray the presence 



of extended caustic structure (closed and sometimes disjoint curves of formally infinite 
magnification on the sky) which are not generated by single point-lenses. 

Binary lenses have already demonstrated their importance to astronomy in several 
ways. Timing the motion of a caustic across a finite background source of known size 
can be used to measure the relative proper motion of the lens and source. This was 
done by several groups who performed a joint analysis ( |Afonso et al. 20U(]| and references 
therein) of their data to demonstrate beyond reasonable doubt that the lensing system of 
event MACHO 98-SMC-l was not a constituent of Galactic dark matter, but a normal 
stellar binary in the Small Magellanic Cloud (SMC) itself. Caustic structure differentially 
magnifies background sources, providing a means of measuring limb-darkening coefficients 
for their stellar atmospheres, as has been done for a K-giant in the Galactic bulge ([Albrow 



|et al. 1999a| ) and an A-dwarf in the SMC ( |Afonso et al. 2000| ). Mathematically, lensing 



planetary systems represent the extreme-mass-ratio subset of binary (or multiple) lenses 
( |Mao fc Paczyhski 1991| ; [Dominik 1999| ). Two groups (PLANET, |Albrow et al. 1998| ; MPS, 



[Rhie et al. 1999a| ) are conducting intense monitoring projects specifically designed to search 
for lensing planets. The efficiency with which planets can be detected must be determined 
separately for each light curve QGaudi fc Sackett 2000|) ; massive planetary companions have 
been ruled out over a wide range of separations for two different microlenses ( |Rhie et"ar 
r999b| ; [Albrow et al. 1999cD . 



The components of all physical binaries rotate about their common center of mass, a 
phenomenon that alters the shape and relative placement of caustic curves on the sky and 
thus, in principle, the form of any light curve generated by a binary microlens ( [Dominik 



19981 ). Nevertheless, rotation has not yet been demonstrated in any binary microlens. 

The unusual light curve of the Galactic bulge microlensing event MACHO 97-BLG-41 
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displays what appear to be source crossings of two disjoint caustic curves separated in time 
by about a month. The event has aheady been discussed by two groups with extensive 
photometric data sets: MACHO/GMAN, who report that they were unable to find a 
static or rotating binary model compatible with their data ( Alcock et al. 1999[ ), and 



MPS/Wise-GMAN ( [Bennett et al. 1999| ), who — based on MPS+MACHO/GMAN data — 
also discount rotation as an explanation for the features in the light curve, concluding that 
a three-body lens (binary + orbiting Jovian planet) is required to produce the light curve 
structure in MACHO 97-BLG-41. Theirs is the first claim in the literature for a planet 
orbiting a normal binary stellar system, as opposed to one component of a binary (|Butler| 
et al. 1997| ; [Cochran et al. 1997|) or a binary pulsar ( [Thorsett et al. 1999| ) . Given the small 



chances of discovering such a system from microlensing observations, Bennett et al. (1999) 
argued that they might well be common. 

Here, we present and analyze the PLANET photometric data set for MACHO 97- 
BLG-41, which is completely independent of the MPS-hMACHO/GMAN data set analyzed 
by Bennett et al. (1999). We also find that no static binary model reproduces our data, 
but unlike Bennett et al. (1999) we conclude that a rotating binary model, with physically 
reahstic parameters, does provide a satisfactory fit to all features in our light curve. 
Since the two data sets are disjoint (and since the MACHO/GMAN data sample key 
caustic-crossing regions that our data do not) we examine whether or not our solution can 
provide a satisfactory fit to their data. We find that our solution (derived without reference 
to any data but our own) is in reasonable agreement with MACHO/GMAN data, so that it 
appears likely that a single rotating-binary solution near the one we have found is consistent 
with all available data. This result obviates the need for a third lensing component (planet) 
to explain the unusual light curve of MACHO 97-BLG-41. 

We describe our data for MACHO 97-BLG-41 and their reduction in §2. Our general 
form for a rotating binary model is given in §3. Our modeling analysis is presented in §4, 
including a demonstration of the non-viability of static binary models, and of the success 
of the rotating binary model in fitting not only our data but also the MACHO/GMAN 
data set for MACHO 97-BLG-41 used in Bennett et al. (1999). In §5 we present the 
physical and kinematical parameters for the binary lens derived from our analysis of our 
data, including the first statistical estimate for the mass and rotational period of a lensing 
binary made possible by our measurement of rotation in this system. We conclude in §6 
that MACHO 97-BLG-41 is the first microlensing event in which binary lens rotation has 
been demonstrated, and that it provides a credible, physically-realistic, and elegant model 
for this complex light curve. 
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2. PLANET Photometric Data for MACHO-97-BLG-041 

An electronic alert of the microlensing event MACHO 97-BLG-41 was issued on 18 
June 1997 by the MACHOQ team. Many such alerts are made each year by the MACHO, 
OGLE0 and EROSQ collaborations, allowing a variety of follow-up observations to be 
undertaken by other teams. The MACHO 97-BLG-41 event was described in the MACHO 
alert as in the direction of the Galactic bulge {a =17:56:20.7, 6 = —28:47:42) and possibly 
of short duration. The PLANET collaboration began monitoring immediately. On 29 June 
1997, MACHO reported that the light curve of the event had flattened, but at a level above 
its baseline value, in a manner consistent with expectations for a binary lens. On 2 July 
1997, PLANET[^ reported an increase of brightness of a few percent, and on 23 July 1997 
predicted that a caustic crossing was likely to occur the following day at JD - 2450000 = 
654 ± 1; this was subsequently confirmed observationally by both MACHO and PLANET. 

The PLANET data set for MACHO 97-BLG-41 consists of photometric measurements 
in the V band from three observing sites and in the / band from five observing sites in 
the Southern Hemisphere during a four- month period directly after the initial alert. In 
addition, six measurements were taken after the event had reached baseline in April 1998 at 
the Canopus 1.0m in Tasmania with the same detector system as used in 1997. The data 
set thus comprises 8 separate light curves. 

After fiat-fielding and bias subtraction, the source star was photometered using the 
image subtraction package ISIS 2.0 ( |Alard 1999| ). In this procedure, pixel values are first 
interpolated to register all frames to an image chosen as the astrometric reference for that 
band and site. A photometric template is then formed for each band and site from several 
good-seeing, low-background frames. ISIS 2.0 uses a non-constant kernel with a fixed 
number of variable parameters to transform the template to each image by performing a 
least squares fit to derive the values of the kernel parameters. Profile-fitting photometry 
is then performed to measure the fiux difference of variable objects between each frame 
and the template by convolving the point spread function (PSF) of the template with 
the best-fitting PSF matching kernel for each of the images. Frame-to-frame differences 
in exposure time, background, transparency, and the size and shape of the point spread 
function are automatically taken into account with this package, but since only flux 



iiT 



MACHO/GMAN Alert Homepage: http : / /darkst ar .astro . Washington . edu 



""^^OGLE Alert Homepage: http://www.astrouw.edu.pl/^ogle/ogle2/ews/ews.html| 



-"^•^EROS Alert Homepage: http://www-dapnia.cea.fr/Spp/Experiences/EROS/alertes.html 
^''PLANET Collaboration Homepage: http://www. astro. rug. nl/^planet 
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differences are measured, the photometric system must be ahgned with a standard system. 
We have done this by performing a hnear regression between ISIS flux differences and 
DoPHOT-reported fluxes (see below), and placing the latter on a standard system using 
standard stars. The ISIS 2.0 package produces a quality parameter to indicate the quality 
of the subtraction, which is a function of the signal-to-noise in the subtracted frame. To 
produce the cleaned data set used for analysis, we require that this quality flag is greater 
than 0.9, eliminating the poorest ~35% of the frames. In addition, we have discarded 9 
other frames for insufficient exposure time, very strong disagreement with data taken at 
nearly the same time on a smoothly varying portion of the light curve, or sensitivity to 
V-band limb-darkening (3 V-band points too near a caustic), which we cannot constrain. 
The result is 325 /-band and 46 V-band measurements, for a grand total of 371 points in 
the high-quality, cleaned ISIS data set. A summary is given in Table 1. 

We have also reduced our data set with the PSF-profile-fitting package DoPHOT 
( pchechter, Mateo, fc Saha 1993|) , using fixed-position catalogs and the average of four or 



more stable, moderately bright and relatively uncrowded stars in the field as a relative fiux 
standard ( |Albrow et al. 199S ). Final cleaned DoPHOT-reduced data sets were formed by 



eliminating all trailed frames, frames that did not process properly, or frames whose image 
quality or signal-to- noise were not adequate to produce reliable (i.e., DoPHOT type 11) 
measurements for the source star or the relative fiux standard. Instrumental magnitudes of 
the SAAO Im DoPHOT system were then calibrated against contemporaneous observations 
of Johnson-Cousins, UBV(RI)c E-region standards ([Menzies et al. 1989|) at the beginning 



of the 1999 observing season. The instrumental magnitudes of other PLANET sites were 
brought to this system in the modeling process by introduction of an alignment parameter 
for every site and band. This parameter is insensitive to widely different models and 
thus has no effect on the results we present here. Linear regression between the ISIS flux 
differences and the calibrated DoPHOT fluxes then allowed calibration of the ISIS dataset. 

The eight ISIS-reduced light curves of MACHO 97-BLG-41 are shown in Figure |I|, in 
which we have overplotted the best-fltting model presented in §4.3 used to photometrically 
align the multi-site data. Due to differences in effective resolution at different sites and 
in different bands, each light curve can be subject to different amounts of non-lensed 
blended light. The model yields a best flt to this blended light so that we can plot the 
magnification of the lensed source light only in Figure |l|. Plotting magnifications in Figure |l| 
also allows presentation of V- and J-band data on the same plot, since the unlensed source 
magnification is by definition equal to unity irrespective of passband.|^ The time scale has 



^^ Technically, the V- and /-band magnification curves can differ in regions such as caustic crossings where 
the passband-dependent source profile is important. Since we have eliminated the few V^-band points in these 
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been expressed in Heliocentric Julian Date (HJD) at the middle of the exposure, with a 
fixed offset: HJD' = HJD— 2450000. We use this timing system throughout. 



2.1. Reddening of the MACHO 97-BLG-41 Field from the CMD 

Figure ^ displays a color-magnitude diagram (CMD) of a 1' x 1' field centered on 
MACHO 97-BLG-41 derived from SAAO Im images. The position of the unmagnified 
source star is shown together with that of the blend (unresolved light along the line of 
sight), as determined from the models presented in § |4.2| . Comparison of this CMD with one 
for a larger 5' x 5' field clearly indicates that the reddening is variable across the field. 

Using our standard-star calibration of this 1' subfield, we find that the center of the red 
clump is at 1^ = 16.20 ± 0.05 and {V — J)ci = 2.55 ± 0.05. Comparing this observed color 
to the dereddened color of (V — /)ci,o = 1.114 derived by Paczyhski et al. (1999), we obtain 
a value of the reddening E{V — I) = 1.44. The reddening law Ai = 1.49-E'(V" — /) of Stanek 
(1996) then yields Aj = 2.15. In §^ we will use this value of reddening and extinction 
together with our best estimate of the color and magnitude of the source from our modeling 
to derive an angular size for the source and thence the relative lens-source proper motion. 



3. Model Parameters for a Rotating Binary Microlensing Event 

A static, binary lensing event of an unblended finite source can be described by the 
seven parameters: d, the angular separation between the two masses as a fraction of the 
angular Einstein radius 9^', q = M2/M1, the mass ratio of the two components; the Einstein 
time ^E = ^'e/z^, where /i is the relative proper motion between the source and the binary 
center of mass; Iq, the time of closest angular approach between the source trajectory and 
the binary center of mass; uq, the angular separation between the source and binary center 
of mass at time to in units of 6e] a, the angle of the source trajectory on the sky relative to 
the binary axis; and p*, the angular size of the source star in units of ^e- 

Throughout, we use the Einstein ring radius of the total mass M, 

2_AGMD^ 



regions, this is not a concern for us here. 



where Dl is the observer-lens distance, D^ is the observer-source distance, and 
Dls = -Ds — -Dl- We center our coordinate system on the center of mass of the binary lens 
such that at time to the lenses are aligned with the x-axis, with the more massive of the 
two. Ml on the right. The angle a is chosen between the positive x-axis and the source 
trajectory so that Mi remains to the right of the source. 

A binary system in which each partner travels on an ellipse about the common center of 
mass is described by the seven parameters of the reduced Kepler problem: the semi-major 
axis a, the period P, the eccentricity e, the phase (at time to, say), and the three angles 
of orientation. However, two of these, a and one of the orientation angles, are already 
contained in the parameterization of a static lens microlensing event (in which distances 
are expressed in 6^ and positions measured relative to the source trajectory), leaving five 
additional parameters specific to binary rotation. If the motion is circular (e = 0), three 
of these canonical parameters remain: P, 0, and the inclination of the orbital plane with 
respect to the sky. 

We model our MACHO 97-BLG-41 data with the simplest rotation model required 
to fit all the observed features, namely rectilinear relative motion of the lens components, 
which represents the next higher order term in lens motion compared to a static binary 
lens. Only the two components of this motion that lie in the sky plane are measurable since 
motion along the line-of-sight does not alter the caustic structure of the binary. We could 
expect to be able to detect the relative acceleration of the binary components only if we 
have a very strong detection of the rectilinear motion and the observations cover a sizable 
fraction of the rotation period. We therefore model the binary rotation by adding only two 
degrees of freedom over that required for the static binary: A^^ and A^, which measure 
the evolution of the separation vector during the course of the observations in the x- and 
y-directions. Under our assumption of rectilinear motion, the separation vector is thus 

d = do + A (^) (2) 

where do = {dQ,0) is the separation vector at time to, and A = (A^^, A^). This formulation 
has the advantage that we need not assume anything about the eccentricity of the orbit. The 
assumption of circular orbits would reduce the complexity of the full orbital solution, but is 
unlikely to be valid for the stellar binaries separated by several AU to which microlensing is 
sensitive ( puquennoy fc Mayor 19?I| ). Regardless of the eccentricity of the orbit, the first 
order term we are attempting to measure is the sky velocity of the separation vector. As we 
explain in § |4.2| , however, it is both inconvenient and time-consuming to parameterize the 
search for rotating solutions in terms of the 9 (7 static -|- 2 rotating) canonical parameters. 



-9- 



We will search for solutions using a different set of parameters, translating our final models 
back to the canonical system at the end of the process. 

In addition to these parameters, for each of the n = 8 light curves, model parameters 
must be assigned for the flux of the source star, Fs and the unlensed background flux, Fb- 
These can vary among sites due to differences in calibration of instrumental magnitudes 
and effective spatial resolution. We also include a parameter rjsee for each light curve to fit 
a linear relation to the systematic trends between measured flux and seeing generally seen 
in our data (cf. [Albrow et al. 1999c|) . The flux of light curve i is thus modeled as 

Fi{t) = Fs,, A{t) + Fb,, + Vsee,i [^FWHM,i(t) - 2"] , (3) 

where ^fwhm.i is the FWHM of the seeing disk, and A is the magnification from the model. 

Finally, since caustic structure can differentially magnify the finite source, we include a 
limb-darkening parameter F/. Our limb-darkening parameter F is related to the canonical 
linear limb-darkening coefficient c through c/ = 3F7/(F/ + 2) (Albrow et al. 1999b, 
Appendix B). Because our V data are not sufficient to constrain Fy, we eliminate this 
parameter from our modeling and the 3 l^-band points that could be affected by it. 

In sum, we require 9 physical parameters for the rotating binary, 3n = 24 for the 
multi-site photometric alignment, blending, and systematic seeing correlations, and 1 for 
source J-band limb darkening, for a grand total of 34 model parameters. We describe how 
we search for solutions in this unwieldy parameter space in §[4.2|. 



4. Light Curve Analysis 

The PLANET light curve of MACHO 97-BLG-41 in Figure | shows two temporally 
sharp features, one at HJD' =619 and a second at HJD' =654. Between these two broad 
features, the magnification is substantially less than A = 3, so that unless the source is 
strongly blended, it cannot be inside a caustic in this region ( jWitt fc Mao 1995] ). 



The leading rise and roll-over at HJD' =650 indicates an approach to one cusp (a sharp 
discontinuous feature in a caustic curve) before the actual crossing near a second cusp of 
the same caustic at HJD' =654. This crossing must have begun after HJD' = 653.21694, 
our last point on the gentle leading shoulder of the second peak. The end of the crossing 
is marked by the ultra-steep decline in brightness near HJD' =654.5. In the interval 
654.50156 <HJD' < 654.53788, the magnification is falling quite rapidly at a linear rate of 
dA/dt = -133.7 ± 1.6 day-^ 
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The first anomaly, which ends near HJD' ~ 619, initially triggered the microlensing 
alert of MACHO 97-BLG-41, then thought to be a normal short duration single-lens event. 
The points in the interval 619.33864 <HJD' < 619.65447 fall smoothly with a slope of 
dA/dt ~ —2.3 day~^ which, although rapid, is considerably slower than the rate seen during 
the caustic crossing near HJD' = 654.5. The concavity and moderate steepness of the light 
curve near HJD'= 619 indicates that the source exits a caustic very close to or over a 
cusp, and approximately parallel to its axis. The crossing must have begun before our first 
SAAO data point at HJD' = 619.34. We note that the concavity of our data in this region 
is inconsistent with the particular three-body solution presented by Bennett et al. (1999), 
though these data alone do not preclude the possibility of a different three-body solution. 



4.1. Failure of Static Models 

The number, shape, and relative orientation of the closed caustic curves generated 
by a binary lens are fixed by its mass ratio q and separation rf, so that information about 
the caustics from the light curve constrain these two quantities (Schneider & Weiss 1986; 
Erdl & Schneider 1993; Dominik 1999). For a given g, wide separations produce two 
diamond (4-cusped) caustics; close separations generate three caustics, one diamond and 
two triangular (3-cusped); intermediate separations produce one 6-cusped caustic. In all 
cases, one of the caustics, called the central caustic, will be near the position of the more 
massive lens. 

Because the second peak in the light curve of MACHO 97-BLG-41 sits atop a region 
of higher background magnification than the first, it is likely to be associated with the 
central caustic. The cusp approach at HJD' =650 between two widely separated crossings 
that are themselves near or over cusps does not support the notion that all three are due 
to the single caustic produced by an intermediate-separation binary. The proximity of the 
cusp-approach shoulder at HJD' =650 to the near-cusp crossing near HJD' =654 constrains 
the size of the central caustic to be rather small compared to the Einstein radius ^e, an 
indication that the mass ratio q is not near unity, d <^ 1, or rf ^ 1. In the last category, 
caustics are symmetric about the line joining them (the binary axis). Since the light curve 
of MACHO 97-BLG-41 is strongly asymmetric about the second caustic crossing, the source 
trajectory must be poorly aligned with the binary axis and thus could not cross the other 
caustic of a wide binary unless the binary were rapidly rotating. On the other hand, the 
triangular caustics of close binaries with mildly unequal masses are naturally positioned at 
an angle with respect to the axis of symmetry of the central caustic. We therefore begin 
by exploring binary models of this type, asking whether the light curve could be produced 
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by the source traversing the magnification pattern of a close binary whose caustics remain 
motionless during the course of the event. 

Strikingly, and despite the fact that our data do not cover either caustic crossing very 
well, we find that no static binary provides a good fit to the data. Static models that 
reproduce the structure of the strongly asymmetric second peak cannot reproduce the first 
caustic region because the source fails to pass through either of the smaller triangular 
caustics at an earlier time. (The topology of the central and triangular caustics can be 
appreciated by examining Fig. |^.) For the same reason, static models that fit the first light 
curve feature are very poor fits to the second. A static binary microlens model simply 
cannot explain the features of our data for MACHO 97-BLG-41. Bennett et al. (1999) came 
to the same conclusion based on their data for this event, which also led them to search for 
a rotating binary solution. Finding no viable rotating model, they concluded that the first 
caustic crossing was not due to one of the binary's own "natural" caustics, but to additional 
caustic structure induced by a planet. We turn now to the viability of a rotating binary 
model for MACHO 97-BLG-41. 



4.2. Success of Rotating Models 

For certain orientations and data samplings, planetary anomalies may be confused with 
those caused by triangular caustics of close binaries with higher mass ratio. Although the 
position of the triangular caustic in static models of MACHO 97-BLG-41 is too misaligned 
with the source trajectory to explain the first caustic peak, its size is compatible with 
the duration of the first peak, in both our data and those used by Bennett et al. (1999). 
We found that models that simply rotate this caustic pattern on the sky yield many 
of the essential features of the light curve. Encouraged by this, and recognizing that 
caustic crossings widely separated in time should offer the best chance for detecting binary 
rotation, we proceeded to a quantitative search for rotating binary models as described in 
the following section. 



4-2.1. Search Strategy 

As we noted in §^, a rotating binary observed for a small fraction of its period can be 
approximated as two masses moving in straight lines about their common center of mass. 
We parameterize this system by beginning with the 7 parameters that characterize a static 
binary event (cio; Q'j^Ej^O; wq, a,p*), and add A^; and Ay to account for the two dimensions 
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of (projected) relative velocity of the binary. Once a model characterized by these 9 
parameters is chosen, the best fit for the additional photometric parameters F^^i, F-Q^i and 
?7see,i can be determined by linear fitting, which requires a negligible amount of computer 
time. Hence, the complexity of the problem is basically determined by the character of the 
9-parameter space. 

In practice, it is quite difficult to search this canonical parameter space by means of 
standard y^ minimization routines, in part because the dimensionality is large, but mainly 
because the x^ surface is highly irregular. A useful approach is to seek an alternative 
empirical parameterization of the problem in which some of the parameters are fixed, or 
at least highly constrained by subsets of the data. For example, Albrow et al. (1999b) 
showed that for static binaries with a well-sampled fold-caustic crossing, it is possible to 
reduce the primary search space from 7 to 5 dimensions by choosing the parameterization 
carefully. Moreover, they were able to restrict the 5-dimensional search to points outside 
of caustics (for which finite source effects can be ignored or handled analytically so that 
the magnification can be evaluated 100 or even 1000 times faster than for caustic-crossing 
points), thereby rendering a 5-dimensional brute- force search tractable. 

Here we develop a similar strategy. Three quantities are relatively well determined by 
cursory inspection of our data: the time the center of the source exits the first caustic tcc,i, 
the time the trailing limb exits the second caustic tcc,2) and the minimum separation from 
the source center relative to the first cusp -Uc,!- We construct a parameterization with these 
three quantities fixed at initial guesses derived from inspection of the data so that the search 
can be restricted to the remaining six parameters. More specifically, by linearly extending 
the light curves from just before and just after the end of the crossing (see e.g. Afonso et al. 
1998), we find that the end of the second crossing is at tec, 2 = ^cc,2 + ^^2 = 654.544 ± 0.001 
days. Here At2 is the radius crossing time of the source over the second caustic so that tcc,2 
is the time the source center egresses. 

We begin our fitting process by assuming that the center of the source passes exactly 
over the first cusp at tcc,i = 619.34, our first data point in the series of SAAO Im points on 
the concave decline. After one iteration, we find that a better estimate is that the center 
passes closest to the cusp of the first caustic at tcc,i = 619.15. The time between caustic 
crossings is then 5t = tcc,2 — tcc,i — 35.18 days, where we have made the evaluation based 
on our final measurement of At2- The six remaining parameters allowed to vary in our 
search are {dmid, Q, h, ^t2, Sd, 69) . Here dmid is the average of the separations di and ^2 at 
the times of the caustic crossings tcc,i and tcc,2, h is the coordinate of the egress point on 
the central diamond caustic relative to the cusp (as measured along the cusp bisector) in 
units of 6e, Sd is the change in the binary separation d during 6t, and 66 is the change in 
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binary orientation during 5t. In our sign convention, 5d> Q and 59 > Q imply a binary lens 
that is separating and rotating counterclockwise on the sky with time. 

We initiate our search for rotating solutions by considering only data from t > 621 
(and encapsulate all the information from the earlier data in the assumption that the source 
crosses the cusp at t = tcc,i). All but 5 of these points are outside the caustics and so require 
only a single call to the lens solver to evaluate the magnification. Hence, the entire light 
curve can be evaluated quickly (~1 second on our system), permitting a very rapid search 
of the 6 dimensional space. For each such model near a minimum in x^, we examine a grid 
of models for the early data points t < 621. This grid is characterized by two parameters, 
namely the two parameters we had held fixed earlier: the time of the first caustic exit and 
the position of the source relative to the cusp at this crossing. A full integration over the 
source is required for most of these early points, but since the search is on a 2-dimensional 
grid, it can be performed in about 4 minutes. For each 6-dimensional model, we add the 
additional x^ from the best 2-dimensional grid to obtain the overall x^ ■ 



4.2.2. Results 

Joint fits to the multi-site multi-band cleaned ISIS data were performed using the 
empirical parameters and discrete-grid strategy described in the previous section. The 
parameters of this minimum were then converted to the canonical parameterization, and 
the minimum checked via a more conventional, but time-consuming, search using a downhill 
simplex method beginning with a seed near the empirical solution. Since the formal errors 
reported by photometric reduction packages often do not represent the true scatter in 
crowded-field data, systematic effects are often present that depend on site and target star. 
We therefore multiply the error bars of each of our 8 light curves by scaling factors cr/aisis 
that force xl/i^i ~ 3) = 1, where Ni is the number of data points in light curve i. (Three 
degrees-of- freedom are already used by the F^^i, F^^i and 77see,i fitting parameters.) These 
error bar scaling factors, given in Table 1, are necessary to ensure that light curves with 
substantially underestimated photometric uncertainties are not given more weight than 
they are due in the joint fit. As an aside, we note that these scaling factors are closer 
to unity for ISIS reductions than for DoPHOT reductions in most cases we have tested, 
apparently indicating that ISIS-reported errors more accurately reflect the true scatter in 
the photometric data. 

Photometric parameters that depend on site, detector, and passband are presented 
in Table 1 for the resulting best-fltting model to the cleaned ISIS data set. In Table 2 
we list the nine empirical geometric fltting parameters and their corresponding canonical 
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counterparts for both data sets. Uncertainties are presented for the fitted parameters; 
canonical parameters are expressed to more digits than are significant since this precision is 
required to reproduce the solution described by the fitted parameter set. 

To specify a solution to a multi-parameter problem to better than Ax^ = 1, one must 
give each parameter i to a precision B^^ , where B = C^^ and Cij is the covariance 
matrix. If the parameters are highly correlated, as they are for the canonical parameters of 
a rotating binary, then B^^ -C Cj/ , so that the individual parameters must be specified 
to much higher precision than the statistical errors. When we convert from the empirical 
to canonical parameters, we find that numerical errors enter that are small compared to 
Cj/ but comparable to B^^ . We therefore minimize x^ for the canonical parameters 
using downhill simplex with the transformed solution as our starting point. This refined 
minimum is given as the last column of Table 2. 

Most of the parameters in the two solutions are very similar, but c/ differs by ~ 0.6a. 
Thus our two minimization routines associated with our two parameterizations agree at 
Ax^ < 1 but not at Ax^ ^ 1, indicating that one or both has numerical errors at this level. 
Since the numerical errors are smaller than the statistical errors, they are of no concern 
here, but would be if the quality of the data could be improved by a factor ten. 

We have computed the full matrix of correlation coefficients for our empirical, geometric 
fitting parameters and have found that the parameters Deltat2 and Tj are anti-correlated 
at the 94% level. This strong anti-correlation arises because the majority of the constraint 
on Ti comes from the 5 SAAO data points near HJD' = 654.5, for which the limb of the 
source was exiting the second caustic. Measurements near the end of a caustic crossing 
are subject to a degeneracy between the source size and limb darkening, such that a large 
source size can be partially compensated for by a smaller amount of limb-darkening, and 
vice versa. This may restrict our ability to predict the light curve shape while the source is 
inside the caustic region, though, as we show in § |4.3| , this restriction is not severe. 



The model light curve produced by our best model is shown in Figure |l| with the 
PLANET data superposed. As the expanded views on the caustic regions show, this 
rotating binary fit is excellent; the solution reproduces our data in every detail. Binary 
lens rotation in MACHO 97-BLG-41 is not only detected, but rather well determined by 
the PLANET data: 6d = -0.070 ± 0.009 and 69 = 5°.61 ± 0°.36 represents ~13% and 
~6% precision (Icr) for the magnitude of the radial and azimuthal motions on the sky, 
respectively. We find no other rotating solutions with similar x^, and suspect that none 
exists since the covariance matrix for the empirical parameters varies smoothly over a very 
wide range in x^ during our search. 
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The residuals of the ISIS cleaned data set from our best model are shown in Figure ^. 
Since we have assumed rectilinear motion of the lens components, systematic residuals with 
time measured from t = (tcc,2 + i^cc,i)/2 might indicate accelerated motion. The absence of 
systematic effects of any kind in the residuals indicates that our assumption of rectilinear 
motion is a good approximation, but also implies that our data for MACHO 97-BLG-41 are 
not sufficient to detect the angular acceleration of the rotating binary components. 

Our best fits indicate that the unblended color and magnitude of the source are 
(y — /)s = 2.51 ± 0.03 and Is = 16.82 ± 0.02, respectively. These uncertainties are 
independent because the first depends only on the ratio of V and / fiux, while the latter 
depends on other model parameters. The blended light, whether from unrelated stars 
along the line of sight or from the lens itself, is quite faint: our best model yields blend 
magnitudes of l^ = 22.8, with errors large enough to make the detection unreliable, and 
Jb = 20.09 ± 0.47. The small amount of blended light places a constraint on the mass of a 
stellar lens: a primary more massive than Mi = 1.4 M© (and hence a binary more massive 
than M = 1.9 M0) at any position between us and the Galactic bulge would produce 
a brighter blend than is observed. A spectrum taken by Lennon et al. (1997) during 
the second crossing may be able to exclude smaller binary masses if the blue (< 4500A) 
unpubhshed portion of the spectrum fails to indicate the fiux expected from less massive 
(hotter) stars. 

Lennon et al. (1997) deduce from their high signal-to- noise spectrum that the source 
star is a cool giant with log g = 3.2 ± 0.3, Tes = 5000 ± 200, and [Fe/H]= -0.2 ± 0.2 dex. 
Interpolating between the model limb-darkening grids of Claret, Diaz-Cordoves, and 
Gimenez (1995), we estimate that such a source star would be expected to have a linear 
limb- darkening coefficient of c/ = 0.56. Limb darkening of the source, which affects the 
SAAO Im data in the sharp decline after the second caustic crossing at HJD' = 654.5 
and (to a lesser extent) the LaSilla data after the first crossing at HJD' = 618.8, can be 
measured from our light curve data alone. We find F/ = 0.42 ± 0.09, corresponding to a 
linear limb-darkening coefficient of c/ = 0.52 ± 0.10, in agreement with atmospheric models 
for a source of the type deduced by Lennon et al. (1997). Our V data do not meaningfully 
constrain Fy. The dereddened source color of {V — /)s,o = 1-07 ± 0.05 we derive in §|| 
together with isochrones of Bertelli et al. (1994) suggests T^s = 4750 ± 150, somewhat 
cooler but consistent with that of Lennon et al. (1997). 
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4.3. Comparison with MACHO/GMAN Data 

Our rotating binary solution makes use of only PLANET data, which unfortunately 
sample only two neighboring points near the first caustic and five points at the very end 
of the second caustic. Nevertheless, since our solution predicts the full light curve over 
these caustic intervals, we compare our prediction with the MACHO/GMAN data for this 
event, kindly provided to us by A. Becker (1999, private communication). These data 
were presented in Alcock et al. (1999) and incorporated into the analysis of Bennett et al. 
(1999), although with two different treatments of the systematic uncertainties. 

We show the comparison between the PLANET model and MACHO/GMAN data for 
the full light curve of the event and for zooms of the first and second caustics in Figure ^ 
In order to determine the magnifications from the flux levels provided to us, we fit the 
fluxes to a function F{t) = A{t)Fs + Fb where Fb is constrained to be non-negative and 
A{t) is the magnification in our model at the time of observation t. The parameters F^ 
and Fb are the minimum required to convert from observed flux to magnification. No 
other fitting parameter was altered to accommodate the MACHO/GMAN data. The light 
curve is calculated using the / band limb-darkening coefficient F/ = 0.42. Although the 
MACHO/GMAN data were taken in other bands, this affects only a few points and those 
by only a small amount. The residuals of the MACHO/GMAN data from our model are 
shown in Figure ^. 

Without accounting for possible systematics in the data through error bar rescaling, 
we find that our model produces a x^ = 2787 for 1557 MACHO/GMAN data points. The 
same model produces a x^ = 2369 for the 1317 MACHO/GMAN data points before or after 
the microlensing event, indicating that the large x^ is likely due to underestimated error 
bars and not due to problems in the model. 

Overall, our model prediction derived from PLANET data is in good agreement with 
the MACHO/GMAN points, even though their data did not enter the fitting process 
in any way, including in the choice of initial model parameters. Even portions of the 
MACHO/GMAN light curve that are devoid of PLANET data are generally well matched 
by the PLANET model. The one exception is the slope of the two Wise Im points at 
HJD' = 654.26, which are separated by 0.004 days. While the model reproduces the 
midpoint of these two points almost exactly, it predicts a change of —0.5%, whereas the 
measured fiux changes by 2.7 ± 0.4%, a formal discrepancy of 8 a, which is discernible but 
not obvious in Figure ^ due to horizontal compression. This confiict cannot be resolved by 
minor adjustments to the model: matching the data would require putting these two points 
in a very different place relative to the caustic structure, which cannot be accomplished 
without significant changes in the model parameters away from the global minimum found 
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using the PLANET data. Hence, if the PLANET model is even approximately correct, the 
errors on these two Wise points must be severely underestimated. 

The MACHO error bars that we have used and shown in Figures ^ and ^ are 
those reported by the SoDoPhot photometry program with an additional 1.4% added in 
quadrature; this is the standard procedure used by the MACHO team to account for 
systematic uncertainties. The uncertainties displayed for the GMAN data from CTIO and 
Wise observatories are those reported from the ALLFRAME photometric package used for 
their reduction, and have not been modified in any way. As we discuss in §^, when fitting 
models to data with mis-estimated errors, rescaling must be done in order to weight one 
data set properly against another to obtain valid statistical error estimates for the derived 
parameters. A. Becker (1999, private communication) reports that in the MACHO/GMAN 
analysis (Alcock et al. 1999) the GMAN error bars were multiplied by a factor of 1.5, while 
D. Bennett (1999, private communication) reports that the MPS analysis (Bennett et al. 
1999) added 1% in quadrature to the uncertainties for GMAN data sets. Since we do not 
have access to the raw data or even to their global characteristics (seeing, background, 
exposure time), we cannot offer an independent judgment. However, because we are not 
fitting models to these data, but only reporting their residuals from our model fit to 
PLANET data, there is no compelling reason to rescale. For completeness, we note that 
the residuals have x^ = 2787 for 1557 data points with MACHO rescaling only (shown), 
and x^ = 2645 and x^ = 2622 under the two GMAN rescaling schemes listed above. Since 
1489 of the total 1557 points are from the MACHO team, the difference in GMAN data 
rescaling has a small effect on the total x^- 

One might ask how the PLANET model is able to fit MACHO/GMAN data in regions 
interior to the caustic where we have no data and where knowledge of the source size 
(relative to the Einstein ring) is crucial. The excellent PLANET coverage over the larger 
anomaly allows determination of q, d^i^, and the trajectory over the second caustic quite 
accurately, so the path of the center of the source through the caustic region is well known. 
Our data at the end of the second caustic crossing fix the time at which the limb of the 
source exits the caustic. Together, the timed trajectories of the source center and its limb 
allow us to determine the source size with a precision sufficient to match MACHO/GMAN 
caustic data. 



5. Physical and Kinematical Parameters of the Lens System 

The proper motion is given by /i = 6^/t^ where 6^, is the angular radius of the source 
and t^ is the time required for the source to travel its own radius. Since the source trajectory 
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and total time for the source to cross the caustic are well- determined by our data, t^, is 
known accurately from our model. We find t* = At2 sin0 = (0.193 ± 0.005) days, where 
At2 is the time for the egressing source radius to cross the second caustic and = 88.2° is 
the angle between the source trajectory and the caustic at this point. 

From our best-fitting solution, we also find the unblended color and magnitude of the 
source, {V — /)s = 2.51 and Is = 16.82, which imply intrinsic colors of {V — /)s,o = 1-07 
and Isfi = 14.67 using the extinction Aj = 2.15 and reddening E{V — I) = 1.44 we derived 
in § p.l| for this field. In order to determine the angular radius of the source, we rederive the 
empirical surface brightness-color relation of van Belle (1999), using the same data set but 
restricted to the domain 2.0 <V — K < 3.5. These data were generously provided to us by 
G. van Belle (1999, private communication) in advance of publication. We obtain 

log(2e,) + V/5 = 1.2885 ± 0.0063 + (0.2226 ± 0.0133) [(K - K)o - 2.823] , (4) 

where the uncertainties are expressed so that they are uncorrelated, and the half-interval 
of the uncertainty corresponds to Ax^ = 1- Note that throughout we have used 6*,, to 
denote angular radius (typical of the microlensing field) not angular diameter (conventional 
in the field of angular-size determinations). We convert our derived color of the source 
{V - /)s,o = 1-07 ± 0.05, to {V - K)s,o = 2.40 ± 0.15 using the calibrations of Bertelh 
et al. (1994) (ATLAS9 models) and Bessel et al. (1998) (ATLAS9 + NMARCS models). 
These calibrations are nearly independent of metallicity and surface gravity. The surface 
brightness-color relation then yields an angular radius of the source of 

^, = (5.56 ±0.54) /ias. (5) 

The error in equation (^) assumes an intrinsic scatter of 8.7% in the van Belle (1999) surface 
brightness-color relation, which we judge from the excess in the scatter over measurement 
uncertainties, and includes the 2.0% standard error about the mean of the restricted 
relationship we have used as well as our measurement uncertainties in J, {V — I), and Aj of 
0.02, 0.03 and 0.08 mag, respectively. Even rather large (~0.5 mag), unrecognized errors in 
the extinction would have a small effect (< 3.5%) on our estimate of 9^ because extinction 
is accompanied by reddening which affects the angular size estimate in the opposite sense. 

We stress that the angular size can be determined much more accurately than the 
physical size (or distance) of the source. We have checked our determination based on van 
Belle's (1999) empirical relation with angular sizes predicted from the isochrones of Bertelli 
et al. (1994). We find that for a range of metallicities (—0.4 <[Fe/H]< 0) and ages (8 Gyr 
< T < 12 Gyr) that reproduce the intrinsic color {V — /)s,o = 1-07 of the source, stellar 
radii and distances (inferred from our I^q = 14.67) vary by almost a factor of two. Despite 
this, the angular radii formed by the ratio of physical radius and distance are extremely 
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stable, within 2% of our slightly modified van Belle (1999) relation. We conclude that 
equation (|) is robust. Our estimate of 6^ is considerably larger than the (2.9 ± 0.7)/ias of 
Bennett et al. (1999); the cause of this is not completely clear. 

Combining our measurements for 6*^, and t* yields 

/i = — = (50±5)kms-ikpc"\ (6) 

which corresponds to (400 ± 40) kms~^ projected onto the Galactic center at 8 kpc. This 
is somewhat higher than would be expected for a bulge source lensed by either a bulge or 
disk lens, although it could be consistent with either. This large proper motion is what 
one would expect for a disk lens lensing a disk source on the far side of the bulge. Since 
the source is close to the Galactic plane (/,6) = (1.32,-1.94), such a geometry is not 
implausible. 



5.1. Mass-Distance Relation 

The Einstein ring radius can be estimated from our model and the proper motion 
estimate above as 

e^ = fitE = (700 ± 70) /ias. (7) 

Hence, from equation (|1|), we obtain the mass-source distance relation 

/ M \ /l-x\/ Ds \-^ , L'l 

t:^ ^^ =0.47 ±0.10, where x = -^. 8 

If the lens were luminous and in the bulge, the total mass would be constrained by 
the mass of the turnoff to be M = (1 + q)Mi ^ 1.5 Mq. Hence, equation (||) implies 
X ^ 0.75. Since lens-source separations are distributed ~ [1 + (1 — x)^/4(/^ + b"^)]^^ where 
2(/^ + 6^)^/^ ~ 0.08 (Gould 2000), this tends to argue against the lens being in the bulge. 
We make this argument more quantitative in the next section. 



5.2. Statistical Estimate of Binary Lens Mass, Distance and Rotational Period 

From the measurements of d^i^, Sd, 69, and 6'e, we can determine the projected 
separation r^ and the projected relative velocity Vp of the binary up to a scale set by the 
lens distance Dl- We use coordinates parallel and transverse to the separation vector dmid 
and consider only leading-order terms in 6d and 66. One then obtains Vp = dmid ^eDi, and 
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Vp = {Sd/St, dmid^O/St) 9eDl- Inserting the parameter values determined from our model, 
we then find 

^'="*"(4l|^)- «.= M-7.6.9)kms-(j^). (9) 

These measurements can be combined with Kepler's Third Law to estimate the mass 
and distance of the lens. To do so, we assume that binaries of the observed mass ratio 
q = 0.34 are distributed along the line of sight in proportion to the density of their 
primaries, that they have a period distribution as measured by Duquennoy & Mayor (1991) 
for G stars in the solar neighborhood, that they are randomly oriented, and that they are 
uniformly distributed in e^, where e is the eccentricity. We assume that the disk population 
has an exponential scale length of 3 kpc, an exponential scale height that varies linearly 
from 300 pc in the solar neighborhood to 100 pc near the Galactic center, and a local 
surface density in stars and brown dwarfs of 35M0pc~^. We further assume that the bulge 
can be approximated as an isothermal sphere with a rotation speed of 220kms~^, truncated 
at 4 kpc from its center. We take the mass function to be dN/dM oc M~^ for both the disk 
and bulge, but with a cutoff of Mi < 1.1 Mq for the bulge. 

We perform a Monte Carlo integration in order to find the relative probability densities 
of the mass, distance, and period of the binary. First, we construct a set of random 
realizations of the orientation, eccentricity, and eccentric anomaly ip (Goldstein 1980) of 
the binary. For each realization, we determine three quantities: Wr = Tp/a, w^ = Vp/{v^), 
and cos$ = r^ • \p/{rpVp), where a is the semi-major axis of the binary and (f^) is its 
time-averaged square velocity. We first ask whether cos$ is near (i.e., within 10%) of 
the observed value of Sd/[(Sdy + {dmidSOYY^'^ = —0.81. If it is not, we skip to the next 
realization. If it is, we use Kepler's Law and the Virial Theorem, GM = a(f^), to find the 
period consistent with the observed parameters 

Next, by eliminating M in the Virial Theorem in favor of 6e, the value of which is 
fixed by equation (0), we obtain 

x2(l -x)= wrw^yl^-^'^ {e^d^,^ [{5df + {d^,d5ef]^ = Q.2QhWrwl (11) 

where we have assumed D^ = 8 kpc in the evaluation. If Wrw"^ > 0.559, this equation has 
no solution and we skip to the next realization. Otherwise, we find the two solutions and 
the corresponding masses from equation (H). If the higher- mass solution has M > 1.9 M© 
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we eliminate it because the lens would have then given rise to detectable blended light 
(which is not observed: see § ^.2| and Table 1). For bulge lenses we eliminate stars above 



the turnoff, M > 1.5Mq (i.e., Mi > I.IMq). Finally, we weight each realization by the 
following factors: 1) dt/dip = 1 — ecos'?/' to take account of the amount of time spent at 
each phase ip, 2) the density of lenses at x, 3) two factors of x, one for the cross section for 
lensing (oc xO-^) and one for the transverse flux of lenses relative to the sources (oc x/x), 4) 
factors to account for the mass function and period distribution as described above, and 
5) factors to account for the transformation from the theoretical quantities (x, M, P) to 
the observables (c^mid , 5d, 66). The most important of these last factors is that for x, which 
scales oc |x(l — x)/{2 — 3x)\, diverging at x ~ 2/3 thus compensating for the small number 
of simulated events near this value. Note that the factors in (3) are actually oc Dj^Oe, but 
since 6-^ is measured and D^ ~ 8 kpc is approximately known, they can be taken as oc x. 

Figure |^ shows the results of this Monte Carlo integration for the lens distance, mass 
and orbital period, respectively. The noise near x = 2/3 and M ~ 0.95 M© is due to the 
low sampling rate and corresponding large transformation factor described above. The 
probability is five times greater that the lens is in the disk than in the bulge. Note that the 
peaks of the probability distributions for the lens mass and distance do not correspond to 
the same solution. If the lens is in the disk, the probability distribution for its total mass 
peaks at that of an M-dwarf binary with M ~ 0.3 Mq, corresponding (from eq. ^) to a lens 
distance of Z^l ~ 3.1 kpc for a source distance of Ds ~ 8 kpc. The peak at Z^l ~ 5.5 kpc of 
the probability distribution for the lens distance, on the other hand, corresponds to a total 
mass M ~ 1.03 M0. The probability distribution for the orbital period is rather sharply 
peaked around 1.5 yr. A firm detection of acceleration in the relative motion of the lens 
components would tighten the constraints on mass, distance and orbital period beyond that 
made possible by our detection of their relative velocity. 

In order for an orbit to be bound, the ratio rv'^/2GM must be less than unity. 
Solutions derived from our Monte Carlo procedure are bound by construction. It can be 
shown that for any such solution, equation (pT]), together with the Virial Theorem and the 
definitions of Wr and Wy, yields a constraint on this ratio for face-on orbits, namely 

r.pvl _x^l-x) 



2GM 0.53 

Since the right-hand side has an upper bound of 0.28 at x = 2/3, all such orbits are quite 
comfortably bound for any value of x. Only orbital phases and inclinations that would 
increase this ratio by a factor of four would result in unbound orbits. 
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6. Conclusions 



We have presented our photometric data set for the complex MACHO 97-BLG-41 
event, consisting of 325 points in the I and 46 points in the V photometric passbands 
taken during the four month period in which the event was in progress. The light curve 
displays two anomalous regions that cannot be fit by a static binary lens, in agreement with 
a similar conclusion by Bennett et al. (1999) modeling an independent data set. Unlike 
Bennett et al. (1999), however, we find that a rotating binary provides an excellent fit to 
our data. We measure changes in the binary separation (in units of the Einstein ring radius) 
of 5d = —0.070 ± 0.009 and in the binary orientation on the sky of 66 = 5°.61 ± 0°.36 during 
the 35.17 days between the separate caustic events, the first clear observational evidence for 
rotation in a binary microlens. 

Our success in finding this rotating solution was facilitated by the high quality of our 
data for this event and by three strategic choices: (1) searching for a solution involving only 
the first-order rectilinear component of the relative motion of the binary on the sky, (2) 
fitting an empirical parameter set that is best constrained by the light curve rather than 
the equivalent, but more non-linear canonical parameter set, and (3) searching on a grid of 
solutions until a gross minimum is found rather than relying on downhill simplex methods 
only to search for approximate minima. 

Bennett et al. (1999) used the MACHO/GMAN data set to reach their conclusion 
that a three-lens system (binary + Jovian planet) was indicated in MACHO 97-BLG-41. In 
contrast, we find that our rotating binary model provides a good fit to the same data with 
fewer additional degrees of freedom, even though the MACHO/GMAN data were not used 
in any way to fix the physical binary parameters and include points in caustic regions that 
our data do not sample well. Furthermore, the particular three-body model presented by 
Bennett et al. (1999) is strongly inconsistent with our data near the first anomaly in the 
light curve of MACHO 97-BLG-41. We conclude that a third body is not required, and 
that binary lens rotation provides a natural, physically-motivated and physically plausible 
explanation of the light curve of MACHO 97-BLG-41. 

Our modeling also produces a limb-darkening parameter for the background K giant 
corresponding to cj = 0.52 ± 0.10, in agreement with models of stellar atmospheres for 
this type of star. Combining our physical and kinematic measurements with reasonable 
assumptions for the distributions in space and separation of Galactic binaries, we are 
able for the first time to derive kinematic probability distributions for the total mass and 
period of a binary microlens; for MACHO 97-BLG-41, these peak at 0.3 M© and 1.5 yr, 
respectively. 
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Table 1. Photometric Parameters fitted to Cleaned PLANET Data for 

MACHO 97-BLG-41 



Data Set 


# Points 


O'/O'ISIS 


1 see 


T b 


T b 


I BAND 












SAAO Im 


97 


1.15 


-0.0270 


16.82 


20.12 


ESO/Dutcli 0.9m 


58 


1.59 


-0.0178 


16.82 


20.33 


Canopus Im 


95 


0.73 


-0.0632 


16.81 


19.44 


CTIO 0.9m 


49 


1.55 


-0.0290 


16.89 


21.57 


Pertfi 0.6m 


26 


1.21 


-0.0141 


16.81 


19.00 


V BAND 












SAAO Im 


14 


1.45 


0.2725 


19.33 


22.84 


ESO/Dutcli 0.9m 


18 


0.92 


-0.1304 


19.38 


22.03 


Canopus Im 


14 


0.34 


-0.1508 


19.43 


19.56 



'^Tlie parameter rj'^^^ is related to the seeing correlation parameter 77 see of S 

by Vsee = Vsee/{Fs + -Fr), and thus has units of arcsec"^. 

^Fluxes Fs and Fb have been converted to standard Johnson V and Cousins 
I magnitudes. 
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Table 2. Geometric Parameters fitted to Cleaned PLANET Data 

for MACHO 97-BLG-41 





Empirical 




CanonicaP 


dmid 


0.5122 (0.0015) 


do 


0.4787 


q 


0.3434 (0.0065) 


q 


0.3434 


tcc,l 


619.18 (0.03) days 


tE 


24.175 days 


tcc,2 


654.544 (0.001) days 


to 


653.4254 days 


Uc,l 


0.0020 (0.0005) 


Mo 


0.07229 


h 


0.028 (0.001) 


a 


112.37 


At2 


0.193 (0.005) days 


P* 


0.00796 


5d 


-0.070 (0.009) 


A. 


-0.0467 


se 


5°.61 (0°.36) 


\ 


0.0366 


T/ 


0.42^(0.09) 


Cl 


0.455^ 



^Precision given in canonical parameters is necessary to 
reproduce the model, but the last two digits are generally uncertain. 

^For the empirical parameters, P/ corresponds to c/ = 0.52±0.10; 
for the canonical parameters, cj corresponds to P/ = 0.358. For a 
discussion, see §4.2.2. 
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Fig. 1. — Our best rotating binary model shown with cleaned PLANET multi-site data in the 
V and / passbands {Green: Dutch/ESO 0.91m, LaSilla, Chile; Red: SAAO Im, Sutherland, 
South Africa; Blue: Canopus Im, Tasmania, Australia; Yellow: CTIO 0.9m, Cerro Tololo, 
Chile; Magenta: Perth 0.6m, Bickley, Western Australia). Rescaled error bars (§ [4.2.2| ) are 
shown, but are generally smaller than the size of the points. Insets are zooms on the two 
caustic crossing regions. Plotted is magnification versus HJD'. 
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Fig. 2. — Calibrated color-magnitude diagram taken from SAAO Im observations of a 1' x 1' 
field centered on MACHO 97-BLG-41. The position of the unmagnified source and much 
fainter blend (as determined from our modeling) are indicated by the large solid dots. The 
dereddening vector connecting the center of the observed red clump to the intrinsic clump 
characterized by Paczyhski et al. (1999) is also shown. 
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Fig. 3. — The caustic topology of our best-fitting rotating model is shown at times near 
the first and second caustic crossings. Units are angular Einstein radii ^e- The straight line 
shows the source trajectory and the arrow its direction of motion. The binary components 
are shown as small dots. As time progresses, the binary rotates counterclockwise (causing the 
caustic pattern to do the same) and the component lenses move closer together (increasing 
the separation between the triangular outer caustics and the central caustic). The zooms 
better illustrate the source (circle) trajectory as it crosses the central caustic (right) and 
triangular caustic (left). The triangular caustic movement during the entire first crossing is 
shown. The source diameter is indicated by the distance between the parallel lines. 
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Fig. 4. — Magnification residuals (data — model) of the cleaned PLANET data set from 
our best-fitting rotating binary model are shown as a function of HJD'. The color coding 
follows that of Fig. |I[ Error bars have been rescaled to account for systematics ( §4.2.2| ). Full 
moon occurred at HJD' = 620 and HJD' = 649. 



-32- 



O 

• I — I 

o 

• 1—1 

• I — I 

a 

cd 



4( 



30 



20 



"1 I r 



"1 I r 



1 \ r 



4 



1 — I — I — I — r 




I I 



40 




616 618 620 622 




Of** 




30 



20 



10 



\ 




\ 653 654 655 



600 



;20 640 

HJD 



660 680 

:50000. 



700 



Fig. 5. — The curves show the magnification of MACHO 97-BLG-41 as a function of 
time derived from the best-fitting model to the PLANET data, and are thus identical 
to those in Fig. 0. The points are data from the MACHO/GMAN collaboration kindly 
provided by A. Becker (1999, private communication). These data points did not enter 
the fit in any way: they are simply superposed on the model curve after undergoing 
a linear transformation to convert from flux to magnification. Color coding indicates 
observing site and passband {Red: MACHO R; Blue: MACHO B; Cyan: CTIO 0.9m R; 
Yellow: MSO 0.8m R; Green: Wise Im R). Insets show the same two zooms displayed in 
Fig. [^. The agreement in regions devoid of PLANET data, near the beginning of the first 
caustic crossing and throughout the bulk of the second crossing, is particularly striking. 
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Fig. 6. — Magnification residuals of tlie MACHO/GMAN data from our best-fitting model 
to the PLANET data only are shown as a function of HJD'. The color coding follows that 
of Fig. ^ MACHO error bars have been multiplied by a customary 1.4 to account for 
systematics. Error bars for GMAN data are those reported by the crowded field reduction 
packages; no attempt has been made here to increase the size of the GMAN error bars 
to account for systematic uncertainties (§ [4.31 ). Full moon occurred at HJD' = 620 and 
HJD' = 649. 
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Fig. 7. — Physical properties of the binary derived from the rotation measurement and other 
observational constraints (see §|5^). Top: Relative probability density of lens distances Dl 
expressed as fraction x of the source distance D^. The upper and lower lines are for disk 
and bulge lenses, respectively. Middle: Same for the total mass of the binary lens. Bottom: 
Same for the orbital period of the binary. To the right of the dashed vertical lines luminous 
(stellar) lenses are excluded, since their mass would then correspond to more blended light 
than is observed. Bulge lenses have a lower distance cutoff (imposed by the simplistic bulge 
model) and an upper mass cutoff (imposed by the mass of bulge turnoff stars) which translate 
into mass and distance limits, respectively, through eq. (H). 



